Dirac-Weyl fermions with arbitrary spin in two-dimensional optical superlattices 
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Dirac-Weyl fermions are massless relativistic particles with a well-defined helicity which arise in the context 
of high-energy physics. Here we propose a quantum simulation of these paradigmatic fermions using multicom- 
ponent ultracold atoms in a two-dimensional square optical lattice. We find that laser-assisted spin-dependent 
hopping, specifically tuned to the (2^+ l)-dimensional representations of the 5u(2) Lie algebra, directly leads to 
a regime where the emerging massless excitations correspond to Dirac-Weyl fermions with arbitrary pseudospin 
s. We show that this platform hosts two different phases: a semimetallic phase that occurs for half-integer s, 
and a metallic phase that contains a flat zero-energy band at integer s. These phases host a variety of interesting 
effects, such as a very rich anomalous quantum Hall effect and a remarkable multirefringent Klein tunneling. 
In addition we show that these effects are directly related to the number of underlying Dirac-Weyl species and 
zero modes. 
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I. INTRODUCTION 

Graphene and topological insulators have recently stimu- 
lated an enormous interest at both the theoretical and exper- 
imental levels (see [1,2] and references therein). These two 
systems share an attractive and unique property, namely, the 
fact that the electronic transport at low energies is not gov- 
erned by the usual Schrodinger equation, but rather by its rel- 
ativistic counterpart: the Dirac equation [3]. In graphene, a 
single layer of carbon atoms densely packed in a honeycomb 
lattice [1], the band structure corresponds to a semimetallic 
phase whose Fermi surface consists of an even number of iso- 
lated points. Interestingly, the low-energy excitations around 
these points display a relativistic dispersion relation, and can 
be thus described by the two-dimensional Dirac Hamiltonian 
for massless fermions. Conversely, in three-dimensional topo- 
logical insulators, which are semiconducting alloys with a 
strong spin-orbit coupling [2], the bulk band structure corre- 
sponds to a gapped insulating phase. Nonetheless, these ma- 
terials also support a robust surface conductivity which can 
be described by an odd number of two-dimensional massless 
Dirac fermions. The recent interest in both graphene and topo- 
logical insulators is two-fold. On the one hand, they provide 
concrete platforms where interesting phenomena were orig- 
inally predicted in a high-energy context, for example, the 
Klein paradox [4], or axion electrodynamics [5]. On the other 
hand, they also foresee novel and useful device applications 
(see for example Ref. [6]). The emergence of massless Dirac 
fermions in solid-state materials is not only an exciting area 
of condensed matter, but is also becoming an exciting topic in 
cold-atom physics (see e.g. [7]). 

In high-energy physics, a relativistic electron is described 
by a spin-^ fermion, an unavoidable fact that fixes the di- 
mensionality of the Dirac spinor [3]. In clear contrast, the 
pseudospin of an emergent relativistic fermion in solid-state 
materials depends on the geometry of the underlying lattice. 
For graphene, the existence of two interpenetrating triangular 
lattices fixes the dimensionality of the Dirac spinor, which in 



this case corresponds to a pseudospin-^ fermion. Recently, 
some researchers have studied alternative lattices, e.g., the ^ 
lattice [8], the line-centered-square (Lieb) lattice [9-11], and 
the Kagome lattice [12], where emergent massless relativistic 
fermions present a pseudospin- 1 structure. These higher-spin 
relativistic fermions show some distinctive features with re- 
spect to graphene 's Dirac fermions, such as all-angle perfect 
tunneling [11], particle localization [10], and the absence of 
the anomalous quantum Hall effect [9]. In view of these re- 
sults, a natural question arises: could we engineer an exper- 
imental setup where massless relativistic particles with arbi- 
trary spin emerge? Moreover, would these higher spins host 
novel effects that have no counterpart in the low- spin cases? 

As we try to increase the pseudospin by modifying the lat- 
tice geometry, the corresponding structures become increas- 
ingly complex. Due to the limited number of two-dimensional 
Bravais lattices and materials that realize such lattices, this 
construction seems doomed to a failure. In this article, we 
propose an alternative approach which is based on the fact that 
these pseudospin structures, which arise from complex lattice 
geometries, can be reproduced by a standard square lattice 
with a matrix hopping [13]. Such a model can be engineered 
with an optical lattice populated by multicomponent ultracold 
atoms with specific state-dependent hoppings. We show that 
such a setup presents a rich playground where the low-energy 
excitations can be described by massless relativistic particles 
with arbitrary spin. Interestingly these fermions are governed 
by the Weyl-like Hamiltonian H = vfS ■ p , which is the mass- 
less version of the Dirac Hamiltonian but can have any spin, 
and whose eigenstates we refer to as Dirac-Weyl fermions in 
this paper. 

More explicitly, we demonstrate that a spin-dependent hop- 
ping, tuned according to the (2^+ 1) -dimensional represen- 
tations of the 5u(2) Lie algebra, directly leads to a regime 
where the low-energy excitations are Dirac-Weyl fermions 
with pseudospin s. We show that this platform hosts two dif- 
ferent phases at half-filling: a semimetallic phase that occurs 
for half-integer s, and a metallic phase that contains a flat zero- 
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energy band at integer s. In the semimetallic phase, we show 
that a Dirac-Weyl fermion with high spin can be described 
as a collection of spin 1/2 counterparts, where each species 
has a different effective speed of light. Accordingly, the 
low-energy transport is characterized by spin 1/2 Dirac-Weyl 
fermions moving at different velocities, an effect known as 
multi-refringence in the field of optics [14]. As a consequence, 
we find exotic tunneling properties across a potential barrier, 
which shows a remarkable multi-refringent Klein paradox. 
Additionally, we study the properties of this exotic Fermi gas 
in the presence of a synthetic magnetic field (see [15] and 
references therein). We find that the Weyl-Landau problem 
can be mapped onto a Dicke Hamiltonian [16], a well-known 
model in quantum optics which describes the interaction of 
an ensemble of two-level atoms with a quantized mode of the 
electromagnetic radiation. From this insightful mapping, we 
derive the exact solution of the Hamiltonian, and predict the 
interesting consequences of an anomalous half-integer quan- 
tum Hall effect [17]. These predictions are confirmed by the 
numerical evaluation of the topological Chern numbers [18] 
and edge-states [19], which directly give the quantum Hall 
sequence. We obtain a general rule describing the quantum 
Hall effect for Dirac-Weyl fermions with arbitrary spin struc- 
tures, generalizing the anomalous quantum Hall effect for 
pseudospin-^ in graphene [1, 19]. 

The paper is organized as follows: in Sec. II, we describe 
the optical lattice setup to simulate the Dirac-Weyl fermions 
with arbitrary spin. In Sec. Ill, we show how to describe a sin- 
gle Dirac-Weyl fermion with high spin as a collection of two- 
dimensional massless Dirac fermions, and also provide topo- 
logical invariants that can be assigned to these excitations. In 
Sec. IV, we present the anomalous quantum Hall effect which 
appears when these exotic particles are subjected to a syn- 
thetic magnetic field. In Sec. V, we explore the transmission 
properties of these excitations in the presence of a potential 
barrier, leading to multirefringent Klein tunneling. In Sec. VI, 
we discuss how to probe various phenomena predicted in the 
main text. We present our conclusions in Sec. VII, and leave 
some technical aspects to the Appendices. In Appendix A, we 
describe in detail how the model Hamiltonian can be simu- 
lated using laser-assisted hopping in optical superlattices. In 
Appendix B, we discuss the analytical solution of the model, 
and the calculation of the topological invariants which charac- 
terize the emerging Dirac-Weyl fermions with arbitrary spin. 
In Appendix C, we present the details on the calculation of the 
Weyl-Landau levels, the zero-energy modes, and their relation 
to a topological index theorem. In Appendix D, we show the 
impact of a flat band on the zero-energy modes. Finally, the 
technical aspects of the multirefringent Klein tunneling are 
presented in Appendix E. 



II. THE MODEL 

Let us consider an ultracold Fermi gas of ^^K atoms trapped 
in the periodic pattern of an optical superlattice [20]. The 
ground state of this atomic gas corresponds to the hyperfine 
manifold with total angular momentum F = 9/2, and cor- 



responding Zeeman sublevels m^r G {— 9/2, ...,9/2}. In this 
work, we shall focus on a subset of N = {2s + 1) internal 
states, which will be referred to as spins, and represented by 
the fermionic creation (annihilation) operators cj^(c^^). Here, 
r stands for the sites of a two-dimensional square superlat- 
tice, and T G {l...A^} is the internal index. As described in 
Appendix A, we assume that the depth of the superlattice is 
so large that the hopping of the atoms due to kinetic energy 
is inhibited. Instead this tunneling shall be laser assisted by 
Raman transitions to certain excited states that are trapped in 
the secondary minima of the superlattice [20]. Accordingly, 
it is possible to control externally the spin-dependent hopping 
of the atoms along the primary minima of the superlattice. In 
the non-interacting limit, which can be accessed by means of 
Feshbach resonances, the Hamiltonian reads 

^ = -LL^[Tv],/,c;:^^^,,c,, + H.c., (1) 

where the hopping amplitude along the V = x^y direction is 
modified by a spin-dependent operator Ty, and we have set 
the lattice spacing to a = 1. In this work, we will study 
a regime where the hopping operators correspond to a N- 
dimensional representation of the 5u(2) Lie algebra, namely, 
Tx = and Ty = Sy, which fulfill the corresponding algebra 
[S,,S±] = ±S± and [S^,S-] = 2S,, where S± = S^±iSy. We 
note that this particular Hamiltonian cannot be obtained from 
an external non-Abelian gauge field [13, 21], which would 
require that T^^^; belong to a Lie group rather than to a Lie 
algebra. In this paper we only consider a two-dimensional op- 
tical lattice, but the setup can also be implemented in a three- 
dimensional one by adding = in the third direction [20] . 



III. DIRAC-WEYL FERMIONS WITH ARBITRARY SPIN 
AND FLAT BANDS 

In this section, we show that the low-energy excitations of 
the 5u(2) Hamiltonian in Eq. (1) correspond to a particular 
type of relativistic particles, the so-csLllcd Dirac-Weyl fermions 
with arbitrary spin [3]. Surprisingly, the spin-dependent hop- 
ping modifies the non-relativistic theory that describes the gas 
of ultracold atoms, and gives rise to emergent Dirac-Weyl 
fermions of arbitrary spin s at low energy. In addition, we 
show that for half-integer spin, each Dirac-Weyl fermion with 
high spin can be decomposed into a collection of Nj:, = s-\-^ 
independent 2+7 spin 1/2 Dirac-Weyl fermions. We note that 
this type of excitations have raised great interest in the con- 
text of graphene [4]. On the other hand, for integer spin, each 
Dirac-Weyl fermion with high spin gives rise to Nj:^ = s spin 
1/2 Dirac-Weyl fermions together with a single zero-energy 
flat band. Recently, the physics of flat bands has also re- 
ceived considerable attention [8, 10, 11]. The atomic system 
discussed here can be used to explore a variety of interesting 
effects at the forefront of condensed-matter and high-energy 
physics. 

Transforming the Hamiltonian in Eq. (1) to momentum 
space, = -jj Y^k^~^^^Ck%, where L is the number of lattice 



sites and k G [—n^n) x [— 7r,;r) lies within the first Brillouin 
zone, one obtains the following A/^-band model 



k 



-^2^v5'vCos(^v), (2) 



where the spinor ^{k) = {ck\^ ...^CkNY contains the fermionic 
operators. Using the properties of the 5u(2) Lie algebra, this 
Hamiltonian can be readily diagonalized (see Appendix B), 
leading to the following spectrum 



Em{k) me{k) =mJ {2txCoskx) 



- [ItyCoskyY, 
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where m = —s^...^s— \^s. We show in Fig. 1 the result- 
ing band structures with A/^ = 2,3,4,5 internal components. 
Further band structures, corresponding to configurations us- 
ing more internal atomic states (N > 5), are easily derived 
from this figure. As can be observed, for N = 2, one recov- 
ers the familiar Dirac cones that arise in graphene [4] . Con- 
versely, for = 3 the Dirac cones are accompanied by a zero- 
energy flat band [8]. When N = 4 and N = 5, we get very 
interesting double-layer cones that host excitations with pseu- 
dospin 3/2 and 2 respectively (cf. below). As shown in Fig. 
1 (c)-(d), these double-layer cones correspond to different ef- 
fective speeds of light, therefore leading to the birefringence 
phenomenon. 

In general, we may conclude that N energy bands touch 
at four highly- symmetric momenta = {^dj^^dy), where 
dy eZ2, which we shall refer to as Dirac points (see Fig. 2). 
At half-filling, the low-energy excitations can be described by 
the following Hamiltonian 

^eff = Y.'¥^p)Hi{p)^^{p), Hiip) = Y^CydySyPy, (4) 

d,p V 

where p = h{k — K^) is the momentum around each Dirac 
point, and we define Cy = 2ty/h. We stress that in the 
isotropic regime Cj = Cy, the effective Hamiltonian for the 
excitations around = ( f ^ f ) corresponds to that of 
the usual Dirac- Weyl fermions. One finds //w = cpA^, 
where we have used the helicity operator As = S • i.e. 
the projection of the spin along the direction of the linear 
momentum [3]. Therefore, the excitations described by 
Eq. (4) can be interpreted as A^v = 4 relativistic Dirac-Weyl 
fermions with arbitrary spin in an underlying anisotropic 
spacetime with an effective speed of light Cj ^ Cy. We stress 
that for N > 2, the Sy matrices do not satisfy a Clifford 
algebra and therefore, the Hamiltonian in Eq. (4) does not 
correspond to a Dirac Hamiltonian. We also note that due 
to the fermion doubling problem [22], the helicity operator 
cannot be globally incorporated into the lattice. 

At this point it is worth emphasizing that in a high-energy 
context, Dirac-Weyl fermions are usually associated to spin- 
1/2 particles. Interestingly, our experimental platform allows 
us to synthesize Dirac-Weyl fermions with any arbitrary 
spin, integer or half-integer. In fact, we find two different 
phases that depend upon the spin s, and can be thus controlled 
experimentally. 



(a) 



E 




(b) 



1 

E 



first Brillouin zone 



first Brillouin zone 



(c) 




first Brillouin zone 




(d) 



E 



first Brillouin zone 

FIG. 1: Energy bands and Dirac-Weyl fermions. Energy spectrum 
E = E{kx,ky) for (a) A^=2 (b) A^=3 (c) A^=4 (d) A^=5. The wave-vector 
k= (kx^ky) belong to the first Brillouin zone. 



a) Semimetallic phases: For semi-integer spin, m / and 
the Fermi surface consists of A^v = 4 isolated points (see 
Fig. 2). It is possible to find an appropriate basis (see Ap- 
pendix B), where each Dirac-Weyl fermion with high spin is 
described by means ofNj)=s-\-^ spin 1/2 counterparts in 2+1 
dimensions. The corresponding Hamiltonian reads 

^eff = L L (cxma^Px^Cyma^Py^ ❖^(/7), (5) 

d,pm>0 

where Cym = 2tym/his the effective speed of light, ay = dyd^ 



4 



are the so-called Dirac matrices in 2+1 expressed in terms of 
the more usual Pauli matrices, and ^m{p) is a two-component 
Dirac- Weyl spinor (see Appendix B for the details). This is 
the Hamiltonian for spin 1/2 massless Dirac fermions . Let 
us stress that in this semimetallic phase, each Dirac-Weyl 
fermion with high spin corresponds to a collection of spin 1/2 
counterparts propagating with a different speed of light (see 
Fig. 1). This will lead to interesting effects, such as the mul- 
tirefringenet Klein tunneling presented in Sec. V below. 

b) Metallic flat-band phases: For integer spin, we find that 
each Dirac-Weyl fermion with high spin describes N^^ = s 
spin 1/2 counterparts in 2+1 dimensions (see Appendix B). 
Besides, in this case, there is always a m = spin compo- 
nent, and thus a zero-energy flat band. We note that the 
presence of this peculiar flat band shall completely modifies 
the properties of the system. Dispersionless energy bands 
have indeed remarkable consequences on single-particle and 
many -body properties [12, 23-26]. In particular, for N = ?> 
our model shares great similarities with the Lieb lattice 
[9, 11], where the low-energy regime is also described by 
a spin-1 Dirac-Weyl Hamiltonian. Note that in this specific 
relativistic band configuration, the fermion doubling imposed 
by the Nielsen-Ninomiya theorem for lattice fermions is 
substituted by the flat band [11, 27]. 

c) Topological invariants: Due to the Hamiltonian particle- 
hole symmetry, the zero eigenvalues are necessarily doubly 
degenerate and could therefore lead to non- trivial Berry 
phases [28]. We have indeed demonstrated, in Appendix B, 
that each of the underlying spin 1/2 Dirac-Weyl fermions 
carries a Berry phase = 27rmsgn{dxdy) , a property also 
familiar from graphene [4]. Note that for half-integer spin, 
the Berry phase is of abnormal parity = ;rmod(27r), re- 
gardless of the Dirac point or Dirac-Weyl fermion. Therefore 
our semimetal phase (s half-integer) has a Berry phase of n 
which is protected by the particle-hole symmetry. 

Furthermore, the particle-hole symmetry also allows for 
a more complete topological description of the semimetallic 
phases. In Appendix B, we show that a particular topological 
invariant can be assigned to each Dirac point. This topological 
index is given by an integer- valued topological charge jl eZ 
introduced in Ref. [29]. In our case, we find that the differ- 
ent Dirac points carry a topological charge with a different 
sign given by i^"^ = Lm>oMm' where {l^ = sgn(44) G Z2. 
Accordingly, we can give a topological characterization of 
each Dirac point, which is /i^ = ±A/d G Z (see Fig. 2). For 
small perturbations that preserve the particle-hole symme- 
try, this semimetallic phase is topologically protected since 
these topological invariants cannot change. On the other hand, 
when the perturbation is strong enough, topological charges of 
opposite sign may annihilate each other, giving rise to inter- 
esting quantum phase transitions that generalize those studied 
in Ref. [30]. In the presence of a weak harmonic confining 
potential, it has been found that the Dirac points and the char- 
acteristic spectrum survives locally in the trap, provided the 
confining potential varies over a length scale much larger than 
the extent of a unit cell [31]. 



IV. WEYL-LANDAU LEVELS AND THE HALF-INTEGER 
QUANTUM HALL EFFECT 

In this Section, we study a peculiar quantum Hall effect 
which occurs when the Dirac-Weyl fermions are subjected 
to an external synthetic magnetic field. The integer quan- 
tum Hall effect is characterized by the perfect quantization 

2 

of the transverse Hall conductivity Oh = v^, where V is an 
integer. This fundamental phenomenon takes place in two- 
dimensional systems subjected to a strong magnetic field and 
exists in the non-interacting limit. Since the ultracold atoms 
are neutral, one needs to mimic the effects of an external 
magnetic field in order to perform a quantum simulation of 
the quantum Hall effect. A simple way to introduce a mag- 
netic field is to rotate the system where the Coriolis force in 
the rotating frame plays the same role as the Lorentz force 
on a charged particle in a uniform magnetic field. Alterna- 
tively, optically induced gauge potentials by means of laser- 
assisted tunneling can be used, see also Appendix A and 
Refs. [15, 20]. These techniques allow us to modify the hop- 
ping of Eq. (1) according to the so-called Peierl's substitu- 
tion Ty TyQ~^^-fv^^^, where A is a synthetic gauge poten- 
tial giving rise to an effective magnetic field B = V xA. As 
shown in recent experiments [32], ultracold atoms can indeed 
be used to investigate the effects of external gauge fields. In 
the following, we investigate the interplay between the Hall 
plateaus and the underlying spin structure. This study offers a 
generalization of the anomalous quantum Hall effect observed 
in graphene, where zero energy modes contribute in a funda- 
mental manner. Here we demonstrate how this peculiar Hall 
sequence is indeed related to the number of Dirac points and 
zero energy modes, for the general case of arbitrary spin struc- 
tures. This study is performed both for the lattice and for the 
continuum limit of our model. 

a) Continuum description: Weyl-Landau levels. For small- 
enough gauge fields, one can show that the Peierl's substi- 
tution leads to the usual minimal coupling performed in the 
effective Hamiltonian (4). The standard procedure is to re- 




FIG. 2: Schematic representation of the Dirac points and their topo- 
logical description. Inside the first Brillouin zone, we find = 4 
isolated points where conical intersections occur. Around these 
points, one finds that the emergent description is that of Dirac- 
Weyl fermions with arbitrary spin. For semi-integer spin s, these 
Dirac-Weyl fermions have an associated Berry phase of abnormal 
parity = 7rmod(27r), and an integer- valued topological charge 
p"^ = ±Nj) e Z, which ensure the robustness of this phase against 
external perturbations. In this figure, we show the double-layered 
cone structure of 5 = 3/2 Dirac-Weyl fermions. 
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FIG. 3: Ladder scheme of the isotropic Weyl-Landau Hamiltonian. 
The system with pseudospin s can be understood as a (2^ + 1) -level 
atom, where transitions between the different spin projections m 
m + 1 (m ^ m — 1) are dressed by the annihilation a (creation a^) of 
a motional quanta. 



place the canonical momentum by a gauge-invariant quan- 
tity p ^ n = p -\- eA, whose compon ents n o longer com- 
mute [nx,!!^;] = —ih^/l^, where /b = \/h/eB is the magnetic 
length. By introducing the bosonic creation- annihilation oper- 
ators = (n^ + in^) where a = {a'^y, the effective Weyl- 
like Hamiltonian in Eq. (4) is recast into 



WL 



■-H^{p^eA)=gd^aS^^gd-aS-^U.c. (6) 



where we have introduced gd± = h{cxdx ± Cydy) / (2a/2/b). 
Interestingly the problem of Dirac-Weyl fermions subjected to 
a magnetic field [33], hereafter referred to as the Weyl-Landau 
levels (WLLs), can be exactly mapped onto the so-called 
Dicke model which is well known from quantum optics [16]. 
This model, which describes the interaction between a collec- 
tion of two-level atoms and a single mode of the quantized 
electromagnetic field, displays a wide range of interesting 
phenomena (see e.g. [34]). In Fig. 3, we represent schemati- 
cally the physical content of the Weyl-Landau Hamiltonian of 
Eq. (6) in the isotropic regime, namely Cj = Cy =\ c, and for 
J = (1, 1). In this situation, = 0, g^^ = hc/Vll^ =: g, 
and the spin raising (lowering) transitions are dressed by the 
annihilation (creation) of motional bosonic quanta. 

In Appendix C, we present a compact way to solve the 
Weyl-Landau Hamiltonian in Eq. (6) for J = (1, 1), in the 
isotropic regime Cj = Cy. This solution is iterative in nature 
and allows us to take the analytical expressions of the WLLs 
from pseudospin s to {s-\- From this method, one can de- 
rive the exact energy spectrum for different pseudospins, such 
as5' = {^,l,|,2} presented in Table I. In Table I, we also find 
that the energy spectrum of the ^ = ^ WLL is analogous to the 
relativistic Landau levels in graphene [1]. We also observe 
the characteristic dependence of the energies, E oc g oc ^/B, 
which is a hallmark that guarantees the relativistic nature of 
the particles. As can be seen in the Table I, this peculiar 
dependence, E oc ^/B, is also fulfilled in higher pseudospin 
cases. In the ^ = 1 case, we observe a couple of particle-hole 
symmetric levels with analogous properties, but also a novel 
zero-energy Landau level which is completely absent in the 
half-integer spin case. The presence of this particular zero- 
energy WLL will have important consequences in the quan- 
tum Hall response of the system. Finally, for ^ = | and s = 2, 
we observe two pairs of particle-hole symmetric levels, which 
is a consequence of the two underlying species of spin 1/2 
Dirac-Weyl fermions assigned to each Dirac point. Note also 
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= ±gVl5^ 


-3Vll n = 3 


En,±,l=±g^J5(2n-3)-3V4n^- 


-\2n+\l 


n>4 


En,±,2 = ^g^^{2n - 3) + 3 V4f72 - 

En,0 = 


-\2n+\l 





TABLE I: Analytical expression of the Weyl-Landau levels. The 
energies corresponding to the Landau levels of Dirac-Weyl fermions 
with different pseudospin ^ = { ^ , 1 , | , 2} are expressed as a function 
of the coupling strength g = hc/\^l^, and the number of motional 
bosonic quanta n. 



that the dependence on the number of motional quanta n gets 
more involved as the pseudospin is increased. 

In addition to the WLLs presented in Table I, we also study 
the presence of certain special topological solutions that oc- 
cur at zero energy (see the details in Appendix C). These so- 
lutions, the so-called zero-energy modes, play a key role in 
the quantum Hall response of the sample and give rise to the 
half-integer anomaly [17]. The underlying topological modes 
contribute with a fractional transverse conductivity (in units 
of e^/h), even in the absence of interactions. In Table II, we 
show that Dirac-Weyl fermions with a half-integer spin s sup- 
port Nz = s-\- ^ topological zero-modes which are protected 
by a topological Atiyah-Singer index theorem [35, 36]. Be- 
sides, they present half the degeneracy of higher Landau lev- 
els, and thus lead to a half-integer anomaly in the quantum 
Hall response. On the other hand, for integer spin s, there 
are also nontopological zero-energy Landau levels that arise 
from the highly degenerate flat band (see Table. I). As argued 
in Appendix C, these zero modes are not related to an index 
theorem and thus, are not protected. This highly-degenerate 
zero-energy band characterizing the integer-spin case is re- 
sponsible for the vanishing of the half-integer anomaly. As 
confirmed numerically in subsection b) (cf. below), the Hall 
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TABLE II: Analytical expression of the zero-energy modes. The 
Weyl-Landau Hamiltonain in Eq. (6) also yields certain zero-energy 
modes whenever the constraints over the number of motional quanta 
presented in Table I are not fulfilled. In this Table, we list these topo- 
logical zero-energy modes for pseudospins s = {^,1,^,2}, where 
\s,m) refers to the spin state, and \n) to the motional Fock state. 

conductivity of the system fulfills 

(^xy = ^^^v + ^ K s half-integer, 

^xy = T^^wV, s integer, 

where V = 0, 1 , . . . determines the different plateaus, and A^v = 
4 is the number of Dirac points. Therefore, as stated above, 
the half-integer anomaly is only valid for the half-integer spin 
Dirac- Weyl fermion. It is also important to note that due to 
the fermion doubling [22], A^v = 4 in our case, the fractional 
character of the Hall sequence is lost. 

b) Lattice description: Computing the Chern numbers. The 
Hall conductivity can be evaluated numerically by diagonal- 
izing the full tight-binding Hamiltonian in Eq. (1) after the 
Peierl's substitution Ty TyC"^^ /v^^^. Considering the stan- 
dard Kubo formula, the Hall conductivity is given by the 
TKNN expression [37] as follows 

= T L ^ch(band), (9) 

occ. bands 

where \ua) are single-particle eigenstates of (1) and where the 
Chern numbers associated to each occupied band, A/ch(band), 
can be efficiently computed using the method of Ref. [18]. 
Here denotes the Fermi energy, which can be tuned in our 
setup by varying the atomic filling. 

Before analyzing the specific Hall plateaus of the systems 
associated to different spin structures, let us draw their energy 
spectra E = E{<$>) as a function of the dimensionless mag- 
netic flux <l> = InBa^. These computed spectra generalize the 
famous Hofstadter butterfly [38]. The fractal butterfly spec- 
tra corresponding to the cases A/^ = 2,3,4,5 are illustrated in 



Fig.4. For N = 2, one recovers the spectrum of the ;r-flux 
model [19]. For N = 3 one observes a spectrum similar to the 
Lieb lattice [9, 39], which highlights the similarity between 
these two models that share a spin-1 configuration. We note 
the existence of a highly-degenerate flat band lying exactly at 
zero energy. For N > 3, the spectra become more complex 
and show complicated overlaps between butterfly-like sub- 
structures. In particular, one identifies two overlapping but- 
terflies in Fig. 4(c)- (d), each of which belongs to one of the 
two species of spin- 1/2 Dirac-Weyl fermions for ^ = 3/2 and 
s = 2. We note that for integer spin (N odd), the central flat 
band at £" = remains robust for all flux <J>. 

The Hall plateaus corresponding to A/^ = 2,3,4,5 are illus- 
trated in Fig. 5. We expect the Hall sequences to be compat- 
ible with the continuum analysis in the low flux regime, and 
we therefore set <I> = 1/51 in this analysis. We also focus on 
the low-energy range, where the description in terms of Weyl- 
Landau levels is valid. First we note that for all the cases (Fig. 
5 (a)-(d)), steps of A/w = 4 are observed in the ranges E^^ <0 
(hole) and E^^ > (particles). The differences between these 
several Hall sequences occur at half-filling (£'f = 0), where 
the flat band and the number of zero modes play a fundamen- 
tal role. For N = 2 (spin-^), one observes a central step of 
Nw^z = 4x1=4, while for A/^ = 4 (spin-|), one gets a cen- 
tral step of Ny^Nz = 4x2 = 8. This numerical results confirm 
the prediction of Eq. (7) based on the analytical expressions 
for the Dirac points, and zero-energy modes derived earlier. 

For general N even (half-integer spin), one indeed obtains 
that the two central plateaus, around half-filling, are given by 
Gjy = ±A/wA/z/2 (in units of e^/h), therefore giving a fun- 
damental signature of the number of Dirac points and zero 
modes. 

For N odd (integer spin), one finds a vanishing contribu- 
tion of the zero modes to the Hall conductivity: a Hall plateau 
corresponding to Gxy=0 is clearly observed in the vicinity of 
= (cf. Fig. 5 (b),(d)). One can understand the vanish- 
ing of the half-integer anomaly as a consequence of the lack 
of an index theorem for the zero modes [40]. Note that this 
general result is in perfect agreement with the Hall sequences 
computed for the ^ [8] and Lieb [9] lattices (i.e. lattices with 
pseudospin 1). 

In order to deepen our understanding of the different Hall 
conductivity plateaus for odd and even N cases, we further in- 
vestigate the associated edge-states, which can be obtained by 
diagonalizing the system in a cylindrical geometry. In other 
words, the topological properties hidden in the bulk (i.e. the 
Chern numbers) may be visible around the boundaries by the 
holographic bulk-to-edge correspondence [19]. This corre- 
spondence is based on the fact that the edge- states carry the 
Hall current along the boundaries of the system. In Fig. 6, we 
show the corresponding energy spectra E = E{k), where ^ is a 
Bloch parameter, for N = 2 (half-integer spin) and N = 3 (in- 
teger spin). Note that Fig. 6 (a) is similar to the spectrum of 
graphene [cf. Fig. 21 in Ref. [1]], and highlights the anoma- 
lous quantum Hall effect that occurs for half-integer spins: 
the contribution of the zero modes at half-filling confirms the 
aforementioned result Gxy{EY^ = 0+) = N^Nz/2 (in units of 
e^/h). In Fig. 6 (a), we observe that each boundary is pop- 
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FIG. 5: Quantum Hall effect in the Dirac-Weyl fermion system. 
The transverse Hall conductivity as a function of the Fermi energy 
Oh{Ey) displays a sequence of plateaus that are associated to an un- 
derlying topological order: (a) = 2 (5 = 1 /2), (b) = 3 (5 = 1), 
= A{s = 3/2) and (d) N=5 {s = 2). Here we set the magnetic 
flux to = 1/51. 



ulated by two edge- states: this is due to the presence of four 
Dirac points in the first Brillouin zone (in contrast with the 
two Dirac points of graphene that lead to a single edge- state). 
Fig. 6 (b) emphasizes the absence of gapless edge-states in 
the first gap above £" = 0, as observed for all the cases corre- 
sponding to odd N. This analysis confirms the general result 
presented in Eq. (7). 

The absence of an anomalous quantum Hall effect for inte- 
ger s is certainly interesting. As discussed above, this effect 
is also manifested by the zero Hall plateau around E^^ = (cf . 
Fig. 5 (b),(d)) or by the absence of visible edge-states stem- 
ming from the zero energy modes (cf. Fig. 6 (b)). The absence 
of an index theorem in this case prevents the robustness of the 
topological zero modes and therefore, it is reasonable to ar- 
gue that they cannot contribute to the Hall conductivity (since 
this physical observable is topologically protected). Further- 
more, the localized properties of the states rooted in the flat 
band [10] could also explain that their associated zero-modes 
would potentially contribute to edge- states with zero velocity 
(i.e. they would not contribute to the Hall conductivity). The 
latter effect is further detailed in Appendix D. 



V. KLEIN MULTI-REFRINGENCE TUNNELLING 



FIG. 4: Hofstadter-like fractal butterflies in the Dirac-Weyl fermion 
system. The energy spectrum E = E{<^>) as a function of the magnetic 
flux shows a fractal butterfly structure: (a) N = 2 (s = j), (b) N = 
3 (s= 1), (c) = 4 (5 = |) and (d) N=5 (s = 2). The paraboUc 
dependence at low magnetic fluxes can be related to the underlying 
relativistic fermions. Note also that for integer spin, one gets a flat 
band exactly located at zero energy. 



It was shown in [4] that due to the coupling of positive 
and negative chirality channels outside and inside a potential 
barrier, quantum tunneling of Dirac particles in graphene be- 
comes highly anisotropic, and the barrier remains perfectly 
transparent for normal incidence. The fermionic and bosonic 
Klein paradox are both discussed in the literature (see [41] 
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k 

FIG. 6: Energy spectrum and current-carrying edge states. Energy 
spectrum E = E{k) for (a) N=2 and (b) A^=3. Here we set <l> = 1/51. 



and references therein). It was shown that while the central 
mechanism for fermionic Klein tunneling is the Pauli exclu- 
sion principle, the key to bosonic Klein tunneling is stimulated 
emission. Jakubsky et al [42] were able to show that for nor- 
mal incidence, the potential can be gauged away by a unitary 
transformation, leaving a free particle dynamics in disguise, 
where a unitary equivalence between the transformed Hamil- 
tonians becomes encoded in a supersymmetry algebra. Evi- 
dence for Klein tunneling has been obtained from graphene 
pn junction [43, 44] and most recently been simulated using 
trapped ions [45]. 

The different helicities carried by the Dirac-Weyl fermions 
can naturally couple inside the barrier, but there is also a pos- 
sibility to transform one helicity to another outside the bar- 
rier, resulting in a remarkable multirefringence phenomenon 
familiar from optics [14]. In Appendix E we present a detailed 
treatment of the Klein tunneling of Dirac-Weyl fermions with 
arbitrary spin. Using the equations introduced there, a numer- 
ical investigation of the double-layered cone structure of spin 
3/2 particles shows that Klein-birefringence is indeed present 
(see Fig. 7 for a schematic overview). The incident energy 
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FIG. 7: Klein birefringent tunneling. Schematic view of a spin 3/2 
Dirac-Weyl fermion incident on a potential barrier (a) while follow- 
ing the outer cone (b). 




FIG. 8: Klein birefringent tunneling. Double transmissions (bire- 
fringence) of spin 3/2 Dirac-Weyl particles from the double-layered 
cone structures as the width and height of the barrier change. In (a), 
the transmission diagram with the increase of width (up two panels) 
and height (down two panels) of potential barrier are shown. Pa- 
rameters used in the simulation: up two panels, 2{Vo —E)/chk = 5 
fixed with the evolution of kD; down two panels with kD = 50 fixed, 
and the evolution of 2{Vo — E)/chk. (b) shows the blowup of initial 
increase with height of the barrier of the down two panels in (a) to 
illustrate the role of evanescent waves, where the cutoff conditions 
for helicities hi = 1/2 and = 3/2 are clearly visable. 



of the particle is E, where the barrier width is D and height 
Vq respectively. The incident particle is chosen to follow the 
outer layer cone. As such, the evanescent waves do not cause 
a cutoff for any helicity outside the barrier, see Fig. 7(b). 

Fig. 8 shows the transmission of spin 3/2 Dirac-Weyl 
fermions when the width and height of the barrier changes. 
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For normal incidence, the barrier is perfectly transparent as 
in graphene [4], see Fig. 8(a). When the incident angle de- 
viates from normal incidence, the component following the 
inner cone shows periodic peaks in the transmission, which 
is the hallmark of birefringence. There are no resonant con- 
ditions, where the barrier would become perfectly transpar- 
ent at certain incident angles, apart from normal incidence. 
For lower barriers, evanescent waves play an important role 
when coupling helicities outside and inside the barrier, as 
shown in the two lowest panels in Fig. 8(a) and Fig. 8(b). 
For a low enough barrier, and for an incident wave beyond 
the critical angle Od = arcsin{{yo — E) / cfikh\)^h\ = 1/2, 
all the helicities coupled inside the barrier are evanescent 
waves, hence there is no transmission for any component as 
shown in Fig. 8(b). For a combination whose critical angle 
is Qci = arcsm{{Vo —E)/cfikh2) and /z2 = 3/2, there is only 
one negative helicity of h\ = \/2 inside the barrier which is 
coupled in the form of a propagating wave. The transmis- 
sion properties are consequently modified dramatically in this 
regime, see Fig. 8(b). These two boundaries stemming from 
the helicities of h\ = 1/2 and /z2 = 3/2 are clearly visible in 
Fig. 8(b). Herein lies the paradox. For low barriers, there is no 
transmission, while for high barriers, the transmission is high. 
For the case of an incident particle following the inner cone, 
there is an additional cutoff condition exerted by evanescent 
waves with helicity /z2 = 3/2 outside the barrier, i.e., beyond 
the critical angle 6c = arcsin{h\/h2) = arcsm{l/3) = 0.34, a 
particle with helicity /z2 = 3/2 cannot be transmitted. 

For n-layered cones we obtain similar results, where inside 
the barrier the presence of evanescent waves will exert a cutoff 
for each helicity, where with increasing barrier height trans- 
mission is gradually allowed for the different helicities. Out- 
side the barrier, the evanescent waves only exert a cutoff when 
a small helicity is transformed to a large one. The transmis- 
sion spectrum typically shows n-refringence, which we refer 
to as Klein multirefringent Tunneling. This multi-refringence 
is a result of the rather unique and non-trivial helicity of these 
Dirac-Weyl fermions. 



VI. DETECTION 

The experimental realisations of the proposed scenarios and 
the detection of the resulting effects are certainly challenging 
but should still be within experimental reach with state of the 
art trapping and manipulation of atomic ensembles. Any de- 
tection scheme capable of resolving the effects discussed in 
this paper will typically have to be able to distinguish between 
spin states, particle density and momentum distribution. We 
briefly outline here some of the possible techniques one can 
use, and refer the reader to [30] for an extensive discussion 
about detection methods of exotic quasi particles in optical 
lattices. 

Regarding the number of Dirac points we note that this can 
be addressed by measuring the atom density close to a zero 
chemical potential [46] . To obtain the number of zero modes 
and the number of Dirac points one can evaluate, or indeed 
measure, the Hall conductivity around half-filling. This could 



be achieved using the Streda formula based on the atomic den- 
sity which can be measured precisely by the in situ individual 
atom detection as in Ref. [47, 48]. In addition, a measurement 
of the atom density as a function of the chemical potential al- 
lows us to map the density of states (DOSs), where the Van 
Hove singularities in the DOSs are directly related to the num- 
ber of layers of the cone structure, see Ref. [13]. By mapping 
out the momentum distribution using atomic angle-resolved 
photoemission spectroscopy (ARPES) Ref. [49], momentum- 
resolved Bragg spectroscopy Ref. [50], or adiabatic release 
Ref. [51], allows us to map the fermi surface and thus the 
location of each Dirac points. 

The existence of a flat band with integer spin can be de- 
tected by its energy dispersion and its related wave function. 
The flat band gives a peak in the DOSs [8], which can be de- 
tected by measuring the atomic density. More importantly, 
the localization properties resulting from the flat band could 
also be detected [10]. This localization can be observed after 
the weak harmonic confinement of the atoms is removed but 
with the optical lattice kept in place: the atoms occupying the 
Dirac cones will fly away fast while the atoms occupying the 
flat band will remain stuck in the immobile flat band states as 
shown in [10]. In our system, this localization property should 
be observed for odd and can be understood by observing the 
vanishing of the Dirac-Weyl spinor components correspond- 
ing to odd cyclotron modes (see for instance Eq.(C5) in Ap- 
pendix C) - in direct analogy with the flat band wave function 
of the Lieb lattice. The vanishing of the components corre- 
sponding to the odd cyclotron modes of the Dirac Weyl spinor 
can also be confirmed by the color resolution strategies sum- 
marized in Ref . [30]. 

The presence of edge states in the bulk gap above half- 
filling Ref. [52-54] is an intriguing concept. For this un- 
doubtedly challenging endeavor one would need to engineer a 
sharp boundary, with a characteristic length on the order of the 
lattice constant, in order to stabilize the presence of topologi- 
cal edge-states within the center of the trap [54] . Loading the 
atoms into the edge states can be achieved with external light 
pulses [53, 54]. In addition, the dynamical structure factor 
5'(q, (o) from light Bragg scattering can also provide a direct 
way to observe the edge states and bulk states as demonstrated 
in Ref. [52]. The lack of edge-states in the bulk gaps around 
half-filling could be an important indicator of the existence of 
a flat band. 

To detect the Klein tunneling the most natural approach 
seems to be designing a potential barrier for the atoms by op- 
tical or magnetic means, and prepare the atoms with a well de- 
fined momentum, then see if particles have tunneled through 
the barrier to the other region. However, a direct confirma- 
tion of the Klein multirefringent tunneling would require a 
launch of a multicomponent mass current. To launch such a 
mass current, several schemes can be used. For example, one 
can connect the optical lattice to two reservoirs with different 
chemical potentials as in [55, 56], exert a static force from a 
tilted optical lattice [57], or by an effective electric field on the 
atoms from the optical dipole force [58]. Measuring the mass 
current across the barrier will consequently reveal the intricate 
dependence on helicity for the tunneling dynamics. 
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VII. CONCLUSIONS AND OUTLOOK 



In this paper, we have proposed the quantum simulation 
of Dirac-Weyl fermions with arbitrary spin by a particular 
laser-assisted tunnehng in optical lattices. By tuning the spin- 
dependent hopping according to the 5u(2) Lie algebra, we can 
assign an arbitrary spin s to these fermions, and go beyond the 
standard spin-^ regime of high-energy physics. 

We have presented a detailed study, both analytically 
and numerically, of several striking aspects of Dirac-Weyl 
fermions. In particular, our system hosts two different phases: 
a semi-metallic phase for half-integer s, and a metallic phase 
that contains a flat zero-energy band for integer spin s. We 
have shown that the low-energy transport in the semi-metallic 
phase is characterized by multirefringent spin 1/2 Dirac-Weyl 
fermions moving at different speeds. As a consequence, 
we also find an exotic Klein tunneling across a potential 
barrier. In the presence of a synthetic magnetic field, we 
have connected the Weyl-Landau problem to the Dicke model 
known from quantum optics. The corresponding Hamilto- 
nian presents a rich structure of Weyl-Landau levels and zero- 
energy modes whose robustness can be related to an index 
theorem for half-integer s, which also includes an anomalous 
half-integer quantum Hall effect. 

Many interesting avenues remain however unexplored. It is 
possible to engineer a mass term by a particular on-site Raman 
transition [13] in order to explore a massive regime and even 
look for a topological insulating phase. Another interesting 
possibility is to engineer a curved spacetime background fol- 
lowing [59], and study the effects of the higher spin of these 
Dirac-Weyl fermions. One can also separate each spin 1/2 
Dirac-Weyl component of the Dirac-Weyl fermions with high 
spin by an appropriate tailoring of the spin-dependent hopping 
[11, 30], which could potentially lead to interesting topologi- 
cal quantum phase transitions. An intriguing supersymmetric 
algebraic structure [42] also deserves further attention in the 
Klein multi-refringent regime. We thus believe that several 
new effects can be explored starting from the results presented 
in this work. 

Acknowledgements. Z. L acknowledges the financial sup- 
port from SUPA (Scottish Universities Physics Alliance). 
A.B. thanks MICINN FIS2009- 10061, CAM QUITEMAD, 
European FET-7 PICC and HIP, UCM-BS GICC-910758, and 
acknowledges very useful discussions with L. Mazza. N.G. 
thanks the F.R.S-F.N.R.S. for financial support, and also V. 
Wens, F. Gerbier, D. Bercioux and J. Dalibard for interest- 
ing discussions. P.O acknowledges support from the Carnegie 
Trust for the Universities of Scotland. We thank R. Shen for 
an interesting discussion regarding the boundary conditions 
for the Klein tunneling. 

Note added. After the submission of this work, another 
manuscript appeared which explored the same idea of Dirac- 
Weyl fermions with high spin but in a different setup [71]. 



Appendix A: Experimental realization of spin-dependent 
hopping 

In this Appendix, we describe the main ingredients of a 
method to engineer the spin-dependent tunneling [20] that 
leads to the Hamiltonian in Eq (1). We consider a cloud of 
ultracold ^^K atoms described by the fermionic field opera- 
tors cl^icjci), where T G {l...A/^} is the internal index that la- 
bels a particular subset of Zeeman sublevels in the ground 
state L = 0,F = 9/2, and x stands for the sites of a two- 
dimensional square superlattice. This particular superlat- 
tice follows from the optical potential created by two pairs 
of counter-propagating lasers along each axis a, V{xa) = 
ViL^ cos^(^L-^a) + V2£^ cos^(2^L-^of), where Vi > V2 repre- 
sent the lattice depths, and an optical wavevector. As 
shown in Fig. 9, the atoms are trapped in a periodic struc- 
ture of primary and secondary minima, that shall be referred 
to as sites r, and links I henceforth. For a deep optical super- 
lattice, the atomic tunneling between neighboring sites will 
be completely suppressed. The fundamental idea to engineer 
a spin-dependent hopping is to assist this tunneling using ad- 
ditional lasers that drive a Raman transition to an excited state 
in the hyperfine manifold F = M = 7 /2 (see Fig. 9), here rep- 
resented by the fermionic operators dl{dx). The use of a pair 
of lasers in a Raman configuration is two-fold. On the one 
hand, the effective frequency can be tuned to the microwave 
transition F = 9/2 ^ F = 7 /2. On the other hand, the effec- 
tive wavevector can be large so that the lasers impart enough 
momentum to the atoms to tunnel between neighboring lattice 
sites. As customary [60], such a two-photon transition is ob- 
tained after the adiabatic elimination of a higher excited state 
in the fine structure L = 1 . The Raman lasers aligned along a 
particular hopping direction, not only drive the transition be- 
tween the hyperfine levels, but also transfer a finite momen- 
tum to the atoms which allows them to tunnel to a neighboring 
site. This assisted hopping [61] is described by the following 
Hamiltonian 

= £ Q.fS„.cl,^d,c-'">^' + H.c, (Al) 

where is the two-photon Rabi frequency driving the tran- 
sition \F = 9/2, t) ^ \F = 7/2, 7/2), and S^^f = (jc|e^^R"|jc') = 
/ d^rw^{r)e^'^^^Wxf{r) determines the momentum transfer that 
the Raman lasers impart on the atom, thus assisting the tran- 
sition between neighboring superlattice sites. The parame- 
ters of this two-photon Raman transition (Or = (Oi — 0)2, and 
kR = ki — k2, follow from each laser frequency and wavevec- 
tor. Here, we have introduced the corresponding Wannier 
functions Wx{r). In the expression of S^^f, one sees the im- 
portance of using a two-photon Raman scheme rather than 
a simple microwave, since the integral between neighboring 
Wannier functions will only be finite when the imparted mo- 
mentum is large (i.e. the effective wavelength is on the order 
of the lattice spacing Ar ^ a, typically a few hundred nanome- 
ters). As shown in [20], by selecting an appropriate detuning, 
Zeeman splitting, and lattice staggering, it is possible to per- 
form an adiabatic elimination of the auxiliary states dl{dx) 
that reside on the links, and thus obtain an effective Hamil- 
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Ionian that describes the hopping of F = 9/2 atoms along the 
primary sites of a two-dimensional lattice. Therefore, the aux- 
iliary link serves as a bus that allows us to assist the tunneling, 
and one obtains the effective Hamiltonian in Eq. (1) 

H = -Y,Y.^v[TyUcl^^^,.Cr, + U.c., (A2) 

where the tunneling strengths ty now depend on the four- 
photon Rabi frequencies. We note that the particular matrices 
Ty can in principle be designed at will, although the experi- 
mental requirements are certainly challenging. 

For the simplest situation s = 1/2, one selects a pair of 
Zeeman sub-levels M = 9/2,7/2, and the hopping opera- 
tors would correspond to = ^Tj^Ty = ^Ty. Each of these 
can be engineered with a single pair of Raman beams, once 
their frequencies are tuned to the resonance given by the Zee- 
man splitting and the lattice staggering. To account for the 
real/complex matrix elements, one should control the laser 
phases appropriately, which can be accomplished by means 
of acusto-optical modulators. The scheme gets more compli- 
cated for ^ = 1, where the hopping operators have now four 
non- vanishing elements, and thus double the number of beams 
required. Let us note, however, that since the different reso- 
nance frequencies rely on the Zeeman splitting and stagger- 
ing, which is a small fraction of the laser frequency [20], the 
desired frequencies can be obtained from the same source, 
once the beams are split, and their frequencies tuned by an 
acusto-optical modulator. For arbitrary s, the scheme is more 
involved, yet benefits from the fact that the hopping matrices 
are of sparse nature. Let us finally remark that there might be 
more clever schemes that take adventage of the light polariza- 
tion to select the different hopping elements. 




L = 



FIG. 9: Superlattice scheme for a laser-assisted hopping. Atoms in 
the ground state L = 0,F = 9/2 are coupled via a two-photon Ra- 
man transition to neighboring atoms in a different hyperfine mani- 
fold L = 0,F = 7/2. This Raman transition takes place through an 
intermediate excited state L= \, such that for large enough detuning 
5, 1 111 1, 1 1 ^ A, where 12/ are the Rabi frequencies for each laser, 
one readily arrives at the Hamiltonian in Eq. (Al). 



Appendix B: Dirac-Weyl fermions and topological invariants 

In this Appendix, we show how to describe a single Dirac- 
Weyl fermion with high spin as a collection of spin 1/2 Dirac- 
Weyl fermions, and also give an explicit derivation of the 
topological charges that can be assigned to these excitations. 
The Hamiltonian in Eq. (2) can be diagonalized using a simi- 
larity transformation for the 5u(2) Lie algebra [62]. In partic- 
ular, the su(2) rotation 

U = S«s^-KSy)^ nl = -2tyC0sky/£{k), (Bl) 

brings the Hamiltonian into a diagonal form 

= U'^HkU = EoPo + Ern{k)G'^. (B2) 

m>0 

where the energies Efn{k) are those of Eq. (3), and we have 
defined Pq = \s,0) {s^0\ sls the projector onto the zero-energy 
flat-band for m = 0, and the Pauli matrices G^ = \s^m){s^m\ — 
\s^ —m){s^ — m|, where \s^m) are the eigenstates of S^. Let us 
note that the eigenstates of the Hamiltonian \Em) = U\s^m) 
can be understood as the 5u(2) spin-coherent states for differ- 
ent fiducial states. 

Once we have an analytic expression for the diagonal 
Hamiltonian, we shall show that for half-integer spin, each 
Dirac-Weyl fermion with high spin can be expressed as a 
collection of spin 1/2 Dirac-Weyl fermions. Let us de- 
fine the Pauli matrices = \s^m){s,—m\ +H.c., and Gm = 
— i|^,m) —m\ +H.C.. Now, we perform the following block- 
diagonal rotation to the Pauli matrices in Eq. (B2) 

U' = diag(e"^?^W(""^^i+"^^i\ ...,e"^?^W(""^^^+"^^^)). 

(B3) 

The block structure of the Hamiltonian (B2), allow us to eas- 
ily derive the effective description in terms of spin 1/2 Dirac- 
Weyl fermions 

Hi/ = U'H^U'^ = - ^ {2tjcmcoskjcG^^2tymcoskyGl^) , 

m>0 

(B4) 

which directly leads to the Weyl-like Hamiltonian in 
Eq. (5). We note that the Dirac-Weyl spinor introduced 
here is built from pairs of states with opposite spin 
^m{p) = U'U(\s^m)^\s^—m)y . By a simple counting 
argument, one finds that there are A^d = + ^ underlying 
spin 1/2 Dirac-Weyl femions per Dirac-Weyl fermion of 
half-integer spin. Conversely, for the integer-spin case, one 
finds N\) = s spin 1/2 Dirac-Weyl fermions and a single 
zero-energy flat band. 

We now derive the Berry phases [28] that can be assigned 
to the semimetallic phase. The parallel transport of the 
eigenstates {\Eyn{k))} below Eq. (B2) can be defined from 
the so-called Berry connection on the principal fiber bundle 
^(T2,t/(1)), where T2 is the 2-torus assigned to the Bril- 
louin zone. The Berry connection is defined as s^m{k) = 
i{Efn{k)\Vk\Em{k)) = i{s^m\ulVkUk\s^m). This expression 
can be evaluated using Sneddon's formula for the parametric 
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differentiation of the operator exponential [62]. In the partic- 
ular case of eigenstates at half-filHng, the Berry connection 
reads s^_\^\{k) = \m\{nlV j^n^ — n^jV kuf) . We can now derive 
the Berry phase 7^ = f^dk- ^rn{k), where C represents a cy- 
cle in the Brillouin zone that surrounds a certain Dirac singu- 
larity. Performing the corresponding line integral, we have 
found that the Berry phase associated to each filled eigen- 
state is = —27r\m\sgn{dxdy) , where dj^dy G {—1, 1} de- 
termine the particular Dirac point. Since the Berry phase is de- 
fined modulo 27r, we find that the Berry phase for each of the 
filled bands has an abnormal parity =71 for half-integer 
spin s. This particular homotopic and topological invariant 
underlies many of the fascinating phenomena of graphene [1]. 

Finally, we close this appendix by calculating an integer- 
value topological charge which is more descriptive than 
the Berry phase presented above. This topological invari- 
ant [29], is defined for particle-hole symmetric Hamiltoni- 
ans, such as Eq. (B4), which fulfills {Hj^^T} = 0, where 
r = diagja^, ...O"^} exploits the block structure of the Hamil- 
tonian and fulfills = F, = I. The topological invariant 
is defined as a winding number /i = ^ /^Tr{F(if^')~Mif^'}, 
which vanishes unless the cycle C encloses a zero of the 
Hamiltonian. Therefore, this invariant is ideally suited to de- 
tect the Dirac points in the Fermi surface, and the value it 
assigns /ij G Z shall be topologically protected for any small 
perturbation that preserves the particle-hole symmetry. This 
winding number reduces to /i^ = £^>o where 

Mm = ^£(dlog/zM-dlog/z;j, (B5) 

and we have introduced h^k = 2txmcoskx-\-2itymcosky. Per- 
forming the line integrals, it can be shown that the topologi- 
cal charges around the Dirac points read fl^ = NY)^gn{dxdy) , 
where A^d = (2^ + l)/2. Accordingly, the charges are integer 
values, and their sign depends on the particular Dirac point 



djc^dy e { — 1,1} under consideration. For particle-hole pre- 
serving perturbations, only when two opposite charges meet, 
a quantum phase transition can take place. Therefore, on very 
general grounds, we can claim that the semimetallic phase is 
topologically protected. 



Appendix C: Analytical solution of the Weyl Landau levels, 
topological zero-energy modes and the Index theorem 

In this Appendix, we derive the exact solution of the Hamil- 
tonian in Eq. (6), and describe the appearance of multiple 
zero-energy modes which are responsible for the half-integer 
anomaly in the quantum Hall sequence. 

Complete analytical solution: In this part of the Appendix, 
we derive the exact solution of the Hamiltonian in Eq. (6) for 
the Dirac point d = (+1,+1), and in the isotropic regime 
Cj = Cy. To do this, we write out the matrix elements ex- 
plicitly as if^L = ^d+E«,m[^^m4,m+lPm^'^ + 

where pn = y/n{2s-\- l—n) and Enm is the (2^ + 1) x (25- + 1) 
matrix with 1 at row n and column m, and zeros otherwise. 
The eigenvector can be expressed as ^ = (0i , 02, • ' ' ? fe+i )^ 
and the eigen-equations are the following operator equa- 
tions pi-ia'^^i-i — e^i + pia^i^i = with / = 1,2, • • • ,2^ + 
1. We solve the Weyl-Landau levels by successive sub- 
stitutions, i.e., we substitute the second equation into the 
first one, and then substitute the third one to the result ob- 
tained by the previous substitution and repeat the process 
(see Table III). After r substitutions, we get Ar^i{aa'^y~^(l)r = 
prBrj{aa'^y~^ a(j)r^i, where Af^i mid Bfj are energy-dependent 
constants, using the Einstein summation convention for re- 
peated indexes ij. Setting pra^(j)r = e0r+i — Pr+i^0r+2 , we 
getAr4a'^ay-\£(t)r^i-pr^ia(t)r^2) =Brjpr{a^ay^r+\' Re- 
casting the expression, we obtain 



[(r+2)/2]>^>/ V [(r+l)/2]>^>i-l V [(r+2) I2]>k>i ^/ 



(CI) 



where (^) is a binomial coefficient and [r] is the integer part of 
r. By using the cutoff condition p2s+\ = , the Weyl-Landau 
levels are obtained compactly as 

[.+3/2] 

A2.+i,K^+l)^"'=0 (C2) 

where A2^+i,i are determined by the recursion above with 
initial values Ai i = e and B\ \ = 1. Here, n is the eigen- 
value of the number operator a^a for the Fock state \n). 
The eigenvector is then expressed as: = {f2s{n,£)\n — 
2s)j2s-i{n,£)\n-2s^l),"'Jo{n,£)\n)y, where fi{n,£) 



are determined by the algebraic equation 

Pi-ifis^i-i^n - 2^ + / - 1 - £f2s^i-i + pifis-iVn -2s ^i = 0, 

(C3) 

with / = 1 , 2, . . . , 2^ + 1 . For ^ = 1/2 and ^ = 1, the results in 
Table III confirm the results in [4, 11]. However, the general 
expression here is applicable to any arbitrary spin. The Weyl- 
Landau states of these high spin particles are a mixture of suc- 
cessive non-relativistic Landau Levels, i.e. the spin and orbital 
degrees of freedom of are highly entangled. This entan- 
glement is a source of many interesting topics, like the meso- 
scopic superposition states in relativistic Landau levels [63], 
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possible fractional Hall states with high Landau-level index 
and non-fractional to fractional QHE transition [64], and non- 
abelian anyons [65]. We note that our Weyl-Landau problem 
can be mapped to a single-mode Dicke model which shows 
a quantum phase transition from the normal to superradiant 
state (i.e. the ground state has a large Weyl-Landau index n) 
in the limit of 5- ^ oo. This may be the way forward to realize 
a stable fractional QHE state with higher Weyl-Landau index 
n in this kind of system when the interaction is turned on. 

Existence of zero-energy modes: In this part of the Ap- 
pendix, we explicitly derive the zero-energy modes displayed 
by the Weyl-Landau Hamiltonian in Eq. (6) for J = (+1, +1) 
in the isotropic regime Cx = Cy. The complete Hilbert space is 
J^ = ^0C^, where ^ = span{|^),^ = 0, 1...} is the Fock 
space of the cyclotron modes, whereas = span{ m) , m = 
—5'..., 5'} is the angular momentum space. Nonetheless, the 
action of the Hamiltonian decomposes in a set of invariant 
subspaces = 0^=o^' where each ^ is spanned by a 



different combination of the Fock and spin states. In particu- 
lar, we find that for n < 2s, there are certain special subspaces 
with odd dimensionality: 

J^o = span{|0)|^,-^)}, 

^2 = span{ |2) 1^, -^) , • • • , |0) + 2)}, 

(C4) 

J^2s-i = span{|2^ - 1)1^, -^), • • • , |0) |^,^ - 1)}. 

Each of these subspaces hosts a zero-energy mode of the 
Weyl-Landau Hamiltonian. By introducing the following 
quantities fn^s.m = 8\/^{s{s-\- 1) —m{m-\- 1)), where g is the 
coupling constant introduced in the Weyl-Landau Hamilto- 
nian of Eq. 6, one finds the following zero-energy modes for 
which if l^o) =0: 



K^) = |o)|.,-.), 

) - |2) 1^, -s) - /2,.,-. |0) 1^, + 2) , 

2s-2 2s-4 
) n fn'^s,s-l-n' -l)\s, s) - fls-l^s-s fl fn',s,s-\-n'\^S - 3) + 2) + • • • + 

w'odd n'odd 

+ (-ir'/Vl,.,.-2 if /2.-l-«V,-.+«'|2)|5,^-3) + (-irl/2 /2,_l_„',,,_,+„'|0)k,5- 1), 

/I'even w'even 

I 



where we have omitted an irrelevant normalization factor. By 
simple counting, we find that the Weyl-Landau Hamiltonian 
hosts Nz = s ^ \ zero-energy modes for half-integer spin. 
Interestingly, the number of zero-modes coincides with the 
number of the spin 1/2 Dirac-Weyl fermions contained by a 
single high spin Dirac-Weyl fermion. Therefore, one may ar- 
gue that each underlying spin 1/2 Dirac-Weyl fermion hosts 
a single zero-energy mode, which shall be responsible for a 
half-integer anomaly in the quantum Hall response of the sys- 
tem (see Eq. 7). Conversely, the number of zero modes for 
integer spin is not bounded (see Table I). 

Index theorem and robustness of zero modes: We have al- 
ready seen two physical manifestations of topology, namely, 
the topological charge, or the abnormal Berry phase, that can 
be assigned to each of the Dirac-Weyl fermions (Sec. Ill), and 
the Chern numbers that determine the quantum Hall response 
of the system (Sec. IV). In this part of the Appendix, we de- 
scribe yet another manifestation of topology: the relation of 
the zero modes to the Atiyah-Singer theorem [35]. This fa- 
mous theorem, which relates the analytical and topological 
features of differential operators, has important consequences 
on the properties of Dirac-Weyl fermions subjected to external 
gauge fields [66]. Graphene therefore has turned out to be an 



excellent platform to understand this relationship both from a 
theoretical [67] and experimental viewpoint [66] . In this Ap- 
pendix, we describe how these concepts can be generalized to 
Dirac-Weyl fermions of arbitrary spin s, and we find that only 
the zero modes of half-integer spin Dirac-Weyl fermions are 
protected by the topological features of the system. As dis- 
cussed in Sec. IV, this justifies the absence of the half-integer 
anomaly for integer- spin Dirac-Weyl fermions. 

The Weyl-Landau Hamiltonian in Eq. (6) for J = (+1, +1) 
in the isotropic regime = Cy presents the following particle- 
hole symmetry 

r, = e''^(^^+^), {H^,r,} = o, (C6) 

which fulfills = I^rJ = r^. TMs operator, known as an in- 
volution [68, 69], allows us to decompose the Hilbert space 
as ^ = J^- , where J^± follow from the orthogonal 
projections P± = ^(I±r^) associated to the ±1 eigenvalues 
of the involution. With this formulation, the Weyl-Landau 
Hamiltonian can be rewritten as a supercharge ^ 
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TABLE III: Weyl-Landau levels by successive substitutions 



n'^ substitution 


expression 


Weyl-Landau levels 


pseudospin 




(^2,1 +A2^2aa'^)(l)2 = P2(52,l)^203 


e = ±gV2n 


5=1/2 


2nd 


(As^l +A3,2^2^^"^)03 = P3(^3,l +B^^2aa^)a^^ ^ 


e = 0,±gV2(2«-l) 


5=1 


yd 


(A4^1 +A4^2^^^^^ +^4,3(<3<3^)^) 04 = P4(54^i -\- B4^2^a'^)^^5 ^ 




1)2 + 9 5 = 3/2 


4th 


(A5,i +A5,2aa"^ +A5,3(aa'^)2) ^5 = (55,1 +Bs^2aa^ +Bs^2>{aa^f) 




-12^+17 5 = 2 





^2,1 = £2 + p2,A2,2 = -p?;52,i = e 

%,i=e(e2+2p2+p2),A3,2 = -e(pi2+p|);53,i=e2+2p2,53,2 = -p? 
^A4,i = + e2(3p2 ^ 2p2 + P32) + ^plp^M.i = -e'(p? + P2 + P3 ) - 
ApIpI,Aa^3 = -p2;54,i =e(e2 + 3p2+2p2),54,2 = -e(pf +p2) 

^A5,i = £(£4 + e2(4p2 + + 2P32 + p2) + 8p2p32 + 4p2p2 + 

3p2p2),A5,2 = e(e2(p2 +p2 +p2 +p2) +6p2p2 ^5p2p2 ^4p2p2)^^^^3 _ 

e(p?p3 +PiP4+P2p4);^5,i =e'(e'+4p? + 3p2+2p32) + 8p2p32,55,2 = 

-e2(p2 +p2 +p2) _ 6p2p32),55,3 = P?p| 



where the differential operators D = P-H^^P^ : Jif^ Jif-, 
and Z)"^ = P^H^^P- : ^ ^ join the orthogonal sub- 
spaces. In this language, the analytical index of the super- 
charge can be expressed as 

ind^ = v+ - v_ = dim(kerZ)) - dim(kerZ)^ ) . (C8) 

For elliptic operators [68], this index can be related to the 
topological features of the system via the famous Atiyah- 
Singer index theorem. This relationship not only gives insight 
into the number of zero modes in the system, but also pin- 
points their robustness with respect to local perturbations of 
the Hamiltonian. From the results in Tables I and II, we ob- 
serve that 

v+ = >y + ^ , v_ = 0, s half-integer, (C9) 

and therefore the total number of zero modes determines the 
index ind=S = s -\- ^, which is related to the total magnetic 
flux that pierces the system. Therefore, these zero modes 
are extremely robust with respect to local perturbations of the 
Hamiltonian. Conversely, the number of zero modes for inte- 
ger spin is unbounded (see Table I). In this case, the differen- 
tial operator D is not an elliptical operator [40], and thus the 
Index theorem does not apply. Accordingly, the zero modes 
for an integer- spin Dirac-Weyl fermion are not topologically 
protected. 

Appendix D: Flat band, zero modes and Hall conductivity 

In Section IV, we have demonstrated that Dirac-Weyl parti- 
cles with integer spin s do not present the anomalous (half- 
integer) quantum Hall effect. Namely, we have shown the 
absence of edge-states stemming from the zero modes, there- 
fore leading to a zero Hall plateau above half-filling. This 
effect is particularly interesting as it seems to be rooted in the 
absence of an index theorem that guaranties the robustness 
of these topological zero modes. An alternative explanation 



comes from the possibility that edge-states with zero veloc- 
ity could be hidden in the flat band and would therefore not 
contribute to the Hall conductivity. 

The squeezing of the edge- states, namely the fact that their 
dispersion relation is given by E{k) =0 for integer s, can 
be investigated in a model that extrapolates continuously be- 
tween the case of Dirac-Weyl particles with integer and half- 
integer spin. Such a model has been introduced by Kennett 
et al. [70]. This model has four sites in its unit cell and con- 
tains a fundamental parameter j3: When j3 = 1, the model is 
equivalent to the Lieb lattice [9], therefore describing a spin-1 
Dirac-Weyl particle and displaying a flat band. For arbitrary 
j3, the energy spectrum displays two interpenetrating cones 
and thus describes a spin-3/2 Dirac-Weyl particle. 

We have computed the Hall conductivity for this model, and 
we find that Gh = ^e^ /h around half-filling for j3 7^ 1. This 
is in perfect agreement with the fact that the model displays a 
single Dirac point. No = I with Nz = 2 zero modes (cf. Eq. 
(7)). For j3 = 1, one observes a zero Hall conductivity plateau 
at half-filling, as expected for a spin-1 Dirac-Weyl particle. 
Therefore, this model is well suited to investigate the fate of 
the edge-states as one goes from Gh = -^e^ /h (half-integer 
spin) to o>/ = (integer spin) around half-filling. 

Thus, we now investigate the edge-states in this model. We 
have computed the energy spectrum of Kennett 's model in the 
presence of an external magnetic field, using a cylindrical ge- 
ometry. For j8 = 0, the first bulk gap hosts a single edge-state. 
As j3 is increased, the energy curves E{k) around E = ±.\ 
are progressively squeezed while opening an energy gap. For 
j8 = 1, this energy gap is exactly located around £ = and 
the energy curves E{k) disappear into the flat band, i.e. they 
are completely squeezed E{k) = 0. Therefore, the edge-states 
associated to these curves, hiding at £^ = 0, have a zero veloc- 
ity in the limit j3 ^ 1, and thus cannot contribute to the Hall 
conductivity. 

This analysis emphasizes the existence of edge-states with 
zero velocity within the flat band, in the presence of a mag- 
netic field. While we have demonstrated this property for 
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Kennett's model, which is equivalent to the Lieb lattice in the 
limit j3 ^ 1 , we beheve that it should be applicable to general 
systems exhibiting Dirac-Weyl particles with integer spin. 

Appendix E: Klein paradox and multi-refringence 

In this Appendix we give a general description of the Klein 
multi-refringence of the Dirac-Weyl fermions with arbitrary 
spin. We shall consider the isotropic, Cj = Cy = c, Weyl-like 
Hamiltonian in Eq. (4), which is valid for low-momentum ex- 
citations around the Dirac point d = (+1, +1), i.e., H^{k) = 
Y^yChSyky. For simplicity, we consider a spin-^ Dirac-Weyl 
particle tunneling through a rectangular potential barrier with 
potential Vq in the interval <x <D and zero elsewhere. The 
particle is incident on the interface at x = at an angle 6 from 
the interface normal (see Fig. 7 ). The plane- wave part of the 
solution is Q^i^xx^kyy)^ where kx = kcos 6 and ky = ksinO. Due 
to the particle-hole symmetry, the helicities form [^+1/2] 
pairs, where [s] stands for the integer part of s. For incident 
particles with energy Vq > E > 0, and helicity ho (i.e. the 
helicity determines the corresponding energy E = cHkho)), all 
negative-helicity channels are coupled inside the potential bar- 
rier with the relation {E — Vq) = cfik^Ji, and k^^ ^ = kf^ cos 
and kf^ ^ = kf^ sin . Outside the potential barrier, helicity ho 
is allowed to convert to other positive values under the condi- 
tion of energy and momentum conservation, i.e., k^^^^h = kho 
and k^^t ,, = ^^^,cos0^^^, and k^^^ y = ^^^,sin0^^,. Due to the 
conservation of parallel wavevectors in the tunneling process, 
one has k^^ sin (j)fj^ = ksinO = ^^^^sin0^^^. Note that a non- 
zero helicity is not allowed to convert into a zero helicity due 
to the violation of energy conservation. The wave function in 
the three regions is: 

h=\orl 

h=—^ or —1 

"Vii = E Ah4e''^''^' + BH<l>t^-'''^''''^\ 0<x<D 

h=—s 
h=^orl 

where the spinor 0^(0^) is the eigenvector of with helic- 
ity h for the right (left) moving wave. Note that and 0^^^^ 
for left moving waves are obtained by using the (;r — 0^^) and 
{n — angles in the solution for the right moving wave. 
Here, r^^A^^Bh and th are 4[^+ 1/2] unknown parameters 
to be determined. By integrating the equation H^^^ = Em 
over an interval in the vicinity of the interface, the bound- 
ary conditions are obtained. For half-integer spin, the bound- 
ary conditions require the continuity of each component of 
the (2^ + l)-component spinor at the two boundaries of x = 



and X = D, which give 4^ + 2 equations that equal to the 
number of unknowns, 4[^+ 1/2] = 4s-\-2 for half-integer 
spin. However, for integer spin where even and odd com- 
ponents of the spinor are decoupled, the boundary condi- 
tions require the continuity of each even spinor component, 
together with the continuity of a sum of two neighbouring 
odd components at x = and x = D (specifically, continu- 
ity of P2i-i^2;-i + p2i^2m with / = 1, . . . ,^). Thesc condi- 
tions give 4^ equations that also equal to the number of un- 
knowns, 4[^ + 1/2] =4^ for integer spin. Thus, in principle, 
the tunneling properties can be completely determined by the 
wavefunction and boundary conditions. Since different he- 
licity spinors carry a different current, the transmission and 
reflection coefficients have to be renormalized with respect 
to the incident current, i.e.,7), = [0^'^^x0^]/[0^^^^x0^^] and 
= thl^L ^x(i>t]/[(i>R^^ Sx(if^^] according to the current density 
j = of the Weyl-like Hamiltonian. Consequently, con- 

servation of the current requires Y!\j2\(^h + ^/i) = 1- In gen- 
eral, for incident particles with spin-^ Dirac-Weyl fermions 
and fixed helicity, the transmission show + 1 /2]-fringence. 

Role of evanescent waves: In graphene, evanescent waves 
don't play a role in the Klein tunneling due to the single- 
layered cone structure. For multiple-layered cones however, 
evanescent waves will play an important role in the coupling 
and transformation of one helicity to another both inside and 
outside the barrier, by exerting cutoff conditions for each he- 
licity component. We first consider the coupling of the inci- 
dent helicity to the ones inside the barrier. For a low enough 
barrier, there are negative helicities inside the barrier which 
are not coupled; they are evanescent. Specifically, for an 
incident wave beyond the critical angle Qc = arcsin{{yo — 
E)/chkm) where m = 1/2 or 1, all helicities are uncoupled, 
thus there is no transmission. For higher barriers, beyond the 
critical angle Qc = arcsin{{Vo —E)/chkhc), only helicities of 
h <hc are coupled in the form of propagating waves. How- 
ever, the transmission properties would be modified dramat- 
ically in this regime. Besides the cutoff condition inside the 
barrier, the evanescent waves also exert cutoff conditions for 
different helicities outside the barrier where a small helicity 
is transformed into a larger one. Since ^s-m^^^^ = sinOh/ho, 
for incident wave beyond the critical angle 6c = arcsm{ho/h), 
the transmitted wave with helicity h become evanescent. This 
gives a cutoff condition for each helicity component which is 
larger than the helicity of the incident particle. It is clear that 
if the incident particle follows the outermost cone of the mul- 
ticone structure, there will be no cutoff for any helicity outside 
the barrier. 

In Section V we present detailed numerical results of Klein 
birefringent tunneling by evaluating the different transmission 
coefficients of spin-3/2 Dirac-Weyl fermions with a double 
cone structure. From the results there, the transmission prop- 
erties of more complex multilayered cone structure can read- 
ily be deduced. 
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